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Abstract. From general considerations of spin-symmetry breaking associated with (anti-)ferromagnetism 
in metallic systems with Coulomb repulsion, we obtain interesting and simple all-order rules involving 
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1 Introduction 

1.1 Theoretical background 

As is obvious and well known [T], magnetic order breaks 
SU(2) spin symmetry (hereafter called SU(2)spin). This 
gives rise to gapless excitations in the form of Nambu- 
Goldstone modes which are magnons and, according to the 
Goldstone theorem, excitations with energy gap, which 
may be called the Higgs bosons. 

These excitations behave like elementary fields, and 
their interaction is central to spin-current conservation 
but, at the same time, they comprise of electrons with 
which they interact: i.e., they arc composite. This imposes 
severe constraints on the properties of these fields, which 
we aim to discuss and exploit in this paper. 

The same situation, of Goldstone fields (i.e., both Gold- 
stone and Higgs fields) that are themselves composite ob- 
jects, arises notably in two problems in the context of high- 
energy physics. The first problem is that of axial symme- 
try breaking at low energy scales due to the SU(3)c strong 
interaction. The Goldstone bosons here are pious. The sec- 
ond problem is that of electro-weak symmetry breaking. 
Although in the Standard Model, the Higgs field, or the 
doublet order-parameter field for SU(2)i(g)U(l)y, is an 
elementary scalar field, it is an attractive possibility that 



this field is composite and, for example, is composed of 
quarks which bind strongly at high scales due to some 
interaction (N.B. not by the SU(3)c' strong interaction 
which is weakly interacting at those scales). The low- 
energy phenomenology would then not be dissimilar to 
that of the Standard Model, but with some constraints on 
quantities such as the Higgs-boson mass. 

A new approach to these problems, due to Gribov, 
have appeared in refs. [51[3J|31[S] . These involve the idea 
of super-criticality and a self-consistent treatment of the 
fermion and Goldstone fields in the presence of the super- 
critical interaction. We shall make use of, and extend, the 
methods presented therein, to the case of magnetism. As 
for the other approaches to these old problems, sec, for 
example, ref. [5] for an old approach to the first problem, 
and ref. [7] for an overview of the various methods and 
techniques developed to handle the second problem. 



1.2 Outline of the paper 

Our work concerns systems of electrons (or holes) which 
interact under a generalized Coulomb exchange (i.e., ex- 
change of a generic gapless photon). We consider the sys- 
tem in the spin-symmetry-broken phase that arise in fer- 
romagnetism and anti-ferromagnetism. 

Our aim is to obtain exact relations between quanti- 
ties that characterize the spin-symmetry-broken phase us- 
ing the Dyson-Schwinger equations. This is possible be- 
cause of the presence of the Ward-Takahashi identities 
which arise because of the conservation of spin symmetry. 
It turns out that the form of the Coulomb interaction does 
not affect these relations. The Coulomb interaction does 
affect, for instance, the electronic self-energy, but these 
arc incorporated in the relations in a general way. 
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Before presenting the full analytical framework, we 
start with the simple case of the discussion of ground- 
state stability in ferromagnetism, in sec. [5] This gives rise 
to an exact rule for ferromagnetism which involves the 
electronic densities of states. We discuss this case with 
illustrations and a physical interpretation. 

The full framework, which employ current-conservation 
techniques, is developed in sees. |3] and HI In sec. El the 
interaction is worked out and presented in the form of 
Feynman rules. In sec. |4l the parameters of the interac- 
tion is worked out. This gives rise to an exact rule for 
anti-ferromagnetism which involves the electronic densi- 
ties of states, but which involves the bare spin exchange 
energy. 

The conclusions are stated at the end. 



2 The ferromagnetic ratio rule 

Before we introduce the full framework, let us discuss the 
stability of the ferromagnetic ground state as an illustra- 
tive example. We do so because this is a relatively simple 
problem, which does not require the full formalism, and 
which nevertheless leads to a strikingly simple and useful 
ratio rule. We shall expand the methods introduced here 
to build the formalism later. 



2.1 Description of the system 

Let us write the electrons in the SU(2)spi,i doublet form: 



(1) 



The Lagrangian has the form: 



er^'Af, ] tpb + (photon K.E.) 



(2) 

Afj^ is the electro-magnetic field, of which we shall retain 
only the electrostatic term Aq later, as the contribution of 
the 3-vector potential is suppressed by the speed of light. 
The photon kinetic energy term may then be taken to be 
— (VAo)^/2, which leads to an electrostatic 1/r interac- 
tion, e refers to the dispersion relation of the electron, e 
is the electro-magnetic charge, which is defined to be neg- 
ative for electrons. F^^ is the vertex function, whose time 
component is 1 for a Lagrangian of this form. 

In the absence of magnetic order, the system is in- 
variant under both the electromagnetic U(1)em and the 
SU(2)spin rotations of ip, where the latter is represented 
by 

U{^,) = exp (iaV,/2) . (3) 

cr' {i = 1,2,3) are the Pauli matrices and 0; are the ro- 
tation angles. Note that we are referring to global rota- 
tions here. Invariance under local rotations requires the 
presence of gauge fields such as the electromagnetic field, 
including the vector potential part and written in a gauge- 
invariant fashion. 



There are two conserved currents, which are orthogo- 
nal. The first is the U(1)em current: 



(4) 



/i refers to the time-space four- vector index (=0,1,2,3). 
The second is the SU(2)spin current. This is written as 



spin 



(5) 



Let us indicate the current diagrammatically by a cross. 
We see that current-current mixing, which we denote as 
^mixing ^'^'^ whose lowcst ordcr term is given by 



(6) 

vanishes by symmetry for all i to all perturbative orders, 
and therefore the two currents are orthogonal. 

When there is magnetic order, there arises, locally, a 
preferred orientation of spin, let us say along J,, and this 
breaks the SU(2)spi„ symmetry, viz: 



SU(2), 



u(i). 



(7) 



As a result, there arises two Goldstone modes whose cou- 
pling is proportional to linear combinations of u^'^, and 
a Higgs mode whose coupling is proportional to u'^ . The 
residual symmetry U(l)z refers to the symmetry under 
rotation by the generator a^: 



Uih) = exp {ia^3/2) = diag (e^^^/^, g-^-^^/^) . (8) 
The form of the effective theory will be discussed later. 

2.2 U(l) current mixing 

A result, which is almost trivial but possibly not pre- 
viously discussed explicitly, is that after this symmetry 
breaking, the currents are no longer orthogonal. Equation 
© is easily calculated. Of particular interest is the — 
component of eqn. ([6|) for i = 3 (i.e., the U(l)z current: 
i = 1,2 vanish) at zero external energy and momenta. The 
vertex function F'^ being equal to 1 when the electrons are 
Fermi-liquid-like for each spin orientation, we obtain 



77°" = lim 

mix A 
g— >0 



Hence 



. f df^ 

im / ; 

^oj (27r 



'^''^^G^{k)G^{q + k)-G^{k)G^{q + k). 

(9) 



77: 



00 



(10) 



Note that 77 is defined with a negative sign, that is 77 = 
—A, where A is the two-point amplitude. G are the elec- 
tron Green's functions. g{tF) are the densities of states at 
Fermi energy. 

As is well known, the Dyson-Schwinger all-order cor- 
rections to a fish diagram such as that indicated in eqn. ^ 
is incorporated by replacing the Green's functions by their 
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all-order counterparts and one or the other of the ver- 
tices (and not both) by their all-order counterpart. Equa- 
tion pH)) is an all-order expression in the sense that the 
Green's functions are arbitrary. However, the vertex cor- 
rection needs care. Let us consider the all-order correc- 
tion to the photon vertex. At zero external energy and 
momenta, the time component of the all-order vertex is 
given by the energy derivative of G~^, by virtue of the 
Ward-Takahashi identity. It follovirs that, if the all-order 
Green's function is given by 



G{E,k) = 



E - e(k) + zOsgn(e(k) - ^i) ' 



(11) 



as is the case for the Fermi liquid, then the vertex is cor- 
rected by On the other hand, there is coming 
from G^, and so the net result is proportional to Z. Z 
being the correct renormalizing factor for the density of 
states, eqn. (jlOl) is exact. It is not difficult to see that a 
more general form of G also admits this property: 



Z(k) 



f{E-e{k)+iOsgIl{e{k)~^l))■' 



(12) 



where / is any function, so long as the density of states is 
definable as the integral of G. Thus it is not a necessary 
condition that the system is a Fermi liquid. 

Returning to eqn. (ITU|) . in general, g-^ and are not 
equal at the Fermi surface, and therefore the currents mix. 

It is worth noting here that the current mixing is zero 
in the case of anti-ferromagnetism, because the two sub- 
lattice contributions are equal and opposite. 

Before proceeding, let us calculate the other fish dia- 
grams. Both for the EM current and for the spin V{l)z 
current, we obtain: 



-'-'em 



n 



00 



(13) 



Again, this is an exact result provided that the Green's 
functions are of the form eqn. (IT2|) and the densities of 
states can be defined as their integrals. Note that although 
we have retained the subscript 'EM' to refer to electromag- 
netism, in reality, we are analyzing electrostatics. 



2.3 Derivation of the ratio rule 

Let us now consider the stability of the ferromagnetic 
ground state. To do so, a primary condition is the van- 
ishing of the tadpole: 




(14) 



as is required by the condition that there are no terms that 
are linear in the Higgs field in the effective Lagrangian. In 
other words, the first derivative of free energy as a function 
of the magnetic order parameter must vanish when the 
ground state is stable. We will also need to check that the 
second derivative is positive. This means the term which is 



bilinear in the Higgs field, or the self-energy of the Higgs 
boson, is positive. The Higgs self-energy is the same as 
77°" calculated earlier on, up to the square of a coupling 
constant. This is necessarily positive. 

Equation is calculated easily, and we obtain 



^tadpole 



(2^) 



.^G,(fc) 



G^(fc)=Pt-Pi- (15) 



p refers to the total density of states of electrons. This 
is an exact expression since higher-order corrections to 
the tadpole are taken into account by making G all-order 
and vertex to be bare. This is always negative if J, is the 
preferred orientation of spin. 

It follows that eqn. (fT4)) by itself is non-zero. How- 
ever, the currents mix, and we must incorporate the con- 
tribution of the photon tadpole, multiplied by the Higgs- 
photon two-point function which has the same form as 



-^mixine Calculated in eqn. pU)) : 



Sab 




(16) 



Now, to make this equation all-order, we must include the 
screening effect in the photon propagator, and this has the 
same form as 77^'^ calculated in eqn. Altogether, we 
obtain 



4 tadpole 
photon part 



= (Pt + Pi)e X 



(.9t(gF) - 9l{<^F))e 



(17) 



The two contributions must vanish when added to- 
gether. Although the charge e appears here, whether one 
takes the charge carriers to be electrons or holes is a mat- 
ter of choice, so e can be taken as constant. Hence, 



gj.(ef) -gt(gF) ^ Pi - Pt 

gd'^p) + gt{<^F) Pi + Pt 



fft('^f ) 



Pi 



(18) 



(19) 



This is our ferromagnetic ratio rule. 

In eqn. ((T8| . the left-hand side is often called the spin 
polarization P, for example in the context of tunnel magneto- 
resistance. The right-hand side is the magnetic moment 
ub — m/fiB divided by the number of carriers n, i.e., 



P = 



(20) 



Note that the definition of n is ambiguous. However, it 
is a measure of the number of electrons or holes that are 
actively involved in the formation of ferromagnetic order. 
As such, one would expect that its order is estimated by 
the number of carriers in the conduction band. If so, we 
obtain a simple rule of the thumb: 



9i > 5t (electrons), 
91 < 9t (holes). 



(21) 
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That is, the density of states of the majority-spin charge 
carriers is always greater. According to this rule, one ex- 
pects the density of states to be rising with energy for 
electrons and decreasing with energy for holes. This can 
be a useful rule of the thumb to establish whether a sys- 
tem with certain given density of states is likely to become 
a ferromagnet. 



2.4 A diagrammatic derivation 

The preceding derivation is formal and, we believe, com- 
plete, but it may appear baffling at the start that the sum 
of apparently only two contributions is sufficient to give 
all-order statements about the system. 

Another way to derive the same result, in a more di- 
agrammatic fashion, is to consider the mixing between 
screened photon and the Higgs boson. Let us denote the 
mixed states by 7 and h. Then the sum of the tadpoles 



7 



■o 



(22) 



needs to be zero, whereas the sum of the tadpoles 

6ab- --4 J 



(23) 



is non-zero, and gives a constant contribution to the en- 
ergy levels of the states. 

The sum over the former set of tadpoles can be ex- 
panded diagrammatically in terms of 7 and h as: 



-O 



■ -cyvx> -o 



(24) 



It will be seen that if the sum of the first two terms is 
zero, then the remaining contributions, which are propor- 
tional to the sum of first two terms, vanish automatically. 
Therefore it suffices to calculate the sum of the first two 
terms. 



2.5 Examples 

One way to understand eqn. ([^ is as a definition of n. 
This number can then be compared with the other esti- 
mates of the number of carriers such as by the Hall effect 
and with the nominal number of electrons or holes. 

As a first example, in the case of half metals, and in the 
ideal case, the spin polarization would be perfect, i.e., P = 
1. The magnetization will also be perfect, and eqn. (j20p 
will be satisfied so long as we take n to be equal to n^. 

Next, in the case of a Coulomb system whose renor- 
malized dispersion relation is given exactly by e(k) = 
(/ik)^/2me, the density of states is given by: 



1 



1 

3^ 



We have defined p as the integral over whole of the oc- 
cupied states. Equation (fT9)) then admits the following 
solutions only: 



(26) 



That is, either the system is a half metal or there is no fer- 
romagnetic order. This is consistent with the observation 
that group 1 elements are not ferromagnetic, and also with 
the observation that the Fermi gas system with quadratic 
(bare) dispersion relation is only weakly ferromagnetic |5] . 

Let us now consider the case of transition metals. In 
eqn. (j20[) . ns is the only quantity which is measured un- 
ambiguously and accurately. P is measurable, but different 
methods yield different results [9,10 . Furthermore, these 
experimental numbers do not match with the results of 
theoretical calculation illi . 

In principle, theoretical numbers for g should be com- 
pared against theoretical numbers for p, and experimental 
numbers for g should be compared against the experimen- 
tal numbers for p. Let us first look at the experimental 
numbers. 



element 


carriers 


Pt 


Pc 


riB 


Fe 


electrons 


0.40 


0.42 ~ 0.46 


2.22 


Co 


holes 


0.35 


0.42 


1.72 


Ni 


holes 


0.23 


0.43 ~ 0.465 


0.606 



element 


n — 


ns/P 


^nominal 


riH 


Fe 


4.8 


~ 5.6 


8 


3.00 


Co 


3.7 


~ 4.9 


3 


0.520 


Ni 


1.3 


~ 2.6 


2 


1.12 



m3e3,/2. (25) 



Table 1. Experimental numbers for spin asymmetry P, mag- 
netic moment ub per site, n calculated as the ratio of these 
two, the nominal number of charge carriers, and the number of 
charge carriers as measured using the Hall effect. Two values 
for P are taken from refs. [9] and [10], respectively, ub is from 
ref. [12]. nH is calculated from Rh listed in ref. [13) . 



In tab. [T] we summarize the experimental numbers for 
P, Ub and n. The calculated values of n are compared 
against the nominal number of charge carriers, i.e., the 
number of 4s and 3d electrons or holes, and with the num- 
ber of carriers calculated from the Hall ratio Rh = —1/ne. 
We see that the the values of n do not seem to be in con- 
tradiction of the nominal number of charge carriers, in the 
sense that n < rinominai for Fe and n w ^nominal for Ni and, 
arguably, Co. However, the variation in the experimental 
numbers is too great to make a concrete statement. 

Let us now turn to the theoretical numbers. 

In tab. [5] we show the numbers for gi/g^ as estimated 
from ref. [IT]. In the case of Fe, p could be estimated 
roughly by the eye as the ratio of the areas underneath 
the density-of-state curves. For Co and Ni, this was not 
possible because of the long tails in these curves. However, 
the relative size of g was consistent with the nature of the 
carriers. That is, the density of states was found to be 
greater for the majority spin. All three cases are thus not 
inconsistent with eqn. (jl9l) . 
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element 






P 


n = ub/P 


Fe 


2.1 ±0.2 


2 ~ 2.5 


0.35 ± 0.04 


6.3 ±0.7 


Co 


7.0 ± 1.5 


> 1 


0.75 ± 0.04 


2.3 ±0.1 


Ni 


10 ± 1.5 


> 1 


0.82 ±0.03 


0.74 ± 0.03 



Table 2. Theoretical numbers for gi/g^ and Pi/ p-^- The num- 
bers for g were measured using a ruler applied to the density- 
of-state plots of ref. [11]. The numbers for p were estimated 
by the eye and a ruler. The Fermi level is at the tail of the 
density of states for the case of Co and Ni, and rendering the 
definition of p difficult or arbitrary. We also show the values of 
P corresponding to these values of g^/g^, and an estimate of 
n based on the actual values of ns which are shown in tab.[T] 



We also show the values of P calculated from g^/g^- 
These are quite different from the experimental numbers 
introduced earlier. As a result, the numbers for n differ 
from before, if we use the same values of n^. 

To summarize, it is difficult to check the ratio rules 
quantitatively, at the present level of accuracy. 



2.6 Physical interpretation 

Let us discuss tadpole cancellation in a more intuitive 
fashion. 




sea ► sea 



Fig. 1. The interaction of a conduction electron with the 
Fermi sea of electrons. 

Figure [T] represents the interaction experienced by a 
conduction electron or holes due to the surrounding Fermi 
sea of (conduction) electrons or holes. Let us say that the 
charge carriers are electrons. Note that fig.H^ corresponds 
to eqn. (ITO . and fig. HJi corresponds to eqn. (ITBl) . 

The interaction shown in Fig.[T]i involves the exchange 
of the density fluctuation of spin, i.e., the Higgs boson. 
This boson being a scalar, its exchange is always repulsive 
between like particles, i.e., between the same spin. Since, 
by definition, there are more majority-spin electrons than 



minority spin electrons, this interaction makes majority- 
spin electrons more energetically unfavourable. That is. 
Fig. [T^, or the density fluctuation of spin, tends to sup- 
press magnetic order. 

The exchange of the Higgs boson is not the only inter- 
action between the conduction electrons and the sea, and 
in Fig. [TJd, we show the Coulomb exchange. This is always 
repulsive, and is of the same magnitude for both type of 
electrons, and so this diagram does not contribute to the 
formation or suppression of magnetic order. 

Figure [T|d by itself is infinite since the photon propaga- 
tor diverges at zero momentum transfer. This is, as usual, 
remedied by the screening effect which is shown in fig. [1]:. 

The screening effect, such as that shown in fig. [T];, 
usually suppresses the charge. This is because a negative 
charge attracts positive charge, and this positive charge 
tends to cancel the negative charge. 

However, the contribution of fig.[IJi requires more thought. 
The sea electrons, which have negative charge, attracts 
positive charge. When this positive charge has the same 
spin as the conduction electron, i.e., when the positive 
charge suppresses the electronic spin which is aligned with 
the spin of the conduction electron, the positive charge at- 
tracts this conduction electron. On the other hand, when 
the positive charge has opposite spin to that of the con- 
duction electron, then the conduction electron is repelled. 

Whether the interaction of fig. [TJi tends to create mag- 
netic order or suppress it depends on which type of spin is 
more likely to be excited, i.e., on the density of states at 
the Fermi surface. This is the meaning of the tadpole can- 
cellation. In other words, the ferromagnetic ground state 
is stable when the interaction due to the fluctuation of 
spin density, which is mediated by the Higgs boson and 
which always suppresses the polarization of spin, is equal 
and opposite to the contribution due to the electrostatic 
polarization of fig-UJi which, depending on circumstances, 
can counteract it. 

2.7 Comparison with the Hubbard model 

When the Coulomb interaction is screened, the interaction 
becomes point-like in the limit of large screening, i.e.. 



9i{e.F) + g\{e-F) 

where C/ is a constant which can be interpreted as the 
on-site Coulomb repulsion U up to a normalization. Let 
us now see what would happen if we were to start from 
a theory which treats the on-site Coulomb repulsion U as 
the starting point, such as the Hubbard model. 

In this case, eqn. (fM)) is unchanged, but eqn. (fT6)) is 
modified to take the following form: 

"•'<^^^^^^ . ,.8) 

Here, as is usual in the Hubbard model, the spin going 
into the fish part must be opposite to the spin going into 
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the tadpole. That is, 



3 Analysis of spin current conservation 



(29) 



in the place of eqn. ((T7)) . This would lead to different con- 
sequences. 

The origin of this discrepancy is clear. Expressed in 
terms of the screened Coulomb propagator, there are two 
contributions that go into eqn. ((28]) . One is the genuine 
tadpole-like contribution of the form eqn. ([T5| . The con- 
tribution of this term is given by 



(30) 



so that this term has the same form as in eqn. (j28p . On 
the other hand, there is a second contribution, which is 
the self-energy correction: 




(31) 



(32) 



The contribution due to this term is given by 

U[gi{eF)pi- gt{<^F)pt] ■ 



Adding together these two contributions yields eqn. (1291) . 

The discrepancy comes because in our approach, the 
self-energy correction is absorbed in the all-order Green's 
function, whereas in the Hubbard-model approach, this 
is not possible. In the Hubbard model, either both con- 
tributions are treated as a tadpole, or both contributions 
are treated as a self-energy correction. If the latter, one 
will have the condition that the simple tadpole, with the 
self-energy corrections, by itself vanishes. This condition 
requires Pf = /Cj,, and therefore we will never have a stable 
ferromagnetic solution out of the Hubbard model. 

This, in our opinion, is a limitation of the Hubbard 
model. The limitation is due to the inability to treat current - 
current mixing, which is the basis of our discussion in this 
section. 

One may still argue that on-site Coulomb repulsion is 
present, physically. In other words, the effective screened 
Coulomb propagator, which ordinarily gives a divergent 
contribution at the origin up to a UV cut-off, is not really 
divergent but only large and finite at the origin. 

If so, this may be thought of as a variant of the UV 
cut-off of the Coulomb propagator -Dphoton, which may be 
parametrized, for example, as 



D 



1 



photon \ 



-k2 - ak4 ' 



(33) 



where a is a parameter (positive or negative). Even if this 
is not permissible as a field theory, it is permissible as 
an UV (Pauli-Villars) regularization procedure. It will be 
seen that such a cut-off does not affect our argument at 
all, since our discussion involves zero momentum photons. 
The electron self-energy will be affected by the UV cut-off, 
but this does not affect our results explicitly. 



Let us now move on to the formalism, which is required if 
we are to go beyond the tadpole-level analysis of the pre- 
ceding section. We adapt Gribov's analysis of axial current 
conservation [3J|3] to the context of spin current conserva- 
tion in systems with partial magnetic order. 

To sum up in one phrase, our goal is to start from 
the Coulombic system, which is defined by eqn. and 
solve it as exactly as possible using the Dyson-Schwinger 
equations, under a number of assumptions. 

The major assumption is that of spontaneous sym- 
metry breaking. If the spin symmetry is broken sponta- 
neously, then the Goldstone theorem guarantees the pres- 
ence of Goldstone and Higgs modes. These modes are, 
in terms of the initial Lagrangian, electronic excitations. 
However, in terms of the effective theory that appears at 
the end, they are elementary, and participate in the con- 
servation of the spin current. This is the main property 
that allows us to solve the Dyson-Schwinger equations. 
The other assumptions, such as the linear or quadratic 
form of the magnon dispersion relation and the constancy 
of exchange energy which are sometimes required, are ap- 
proximations, which we believe are viable, that can be 
lifted if one has the computational resources. 

The resulting effective Lagrangian is found to have the 
following interaction term: 



(34) 



Here is essentially the order-parameter field, but with a 
certain formal difference which we shall discuss. We would 
like to emphasize at this point that this equation is not 
our starting point. It is rather the end product of solving 
the Coulombic system by means of the Dyson-Schwinger 
equations, with the aid of the Goldstone theorem and the 
Ward-Takahashi identities. 

Let us start by discussing current conservation. As dis- 
cussed in the previous section, current is absolutely con- 
served when the spin symmetry is conserved. 



Q^f, spin 



0. 



(35) 



Current conservation is refiected in the following Ward- 
Takahashi identity: 



(36) 



is the vertex in the momentum space. Ai^2 refer to 
the spin states, but these are dummy indices here in the 
sense that G^^{q) is independent of A. Thus eqn. (1551) 
holds for any combination of spin, and therefore current 
is conserved, q are d + 1-vectors with components (go, q). 
qo is the energy and q is the spatial momentum, with 

The Ward-Takahashi identity is violated in the symmetry- 
broken phase, since there is now an energy difference AE, 
which is the exchange energy, between the different spin 
states: 

AE = Gl\q)-G:^\q). (37) 



Kosuke Odagiri; Current conservation and ratio rules in magnetic metals with Coulomb repulsion 



7 



This is the definition for the ferromagnetic case, when AE 
is positive if we take J, to be the majority spin state. In the 
anti-ferromagnetic case, AE is positive in one sub-lattice 
and negative in the other. 

If G^^ and G^^ are both linear in energy, AE is given 
by e-f- — and is constant up to a possible dependence 
on the spatial momentum q. In principle, AE depends 
on go and q. In particular, at the threshold, G~^ would, 
in general, have a singular structure corresponding to the 
emission and absorption of the Goldstone boson (jj through 
the process e* — ^ ecfi. The results of the previous section 
are stable against such corrections, as we have discussed. 
However, the results of this section are more easily derived 
for constant AE, which corresponds to the case of the 
Fermi liquid whose exchange energy is constant. 

After the symmetry violation, the currents Jf 2 are no 
longer conserved, and the Ward-Takahashi identity is vi- 
olated by 

r^^{qi-q2)^ « G^^{qi)~G^^{q2)±AE (Ai ^ A2). (38) 

This AE contribution is of the same form as the cou- 
pling of the Goldstone boson (j>, and current conservation 
is restored by including the contribution of the Goldstone 
boson. This is the case even when AE is not constant. 
Specifically, spin current conservation is restored for the 
vertex F which is modified by the inclusion of the Gold- 
stone boson, 



X 
I 



(39) 

As before, the crosses indicate the spin-current vertices, 
and the dashed line indicates the Goldstone boson. There 
is nothing strange in this result, since the Goldstone boson 
arose in the first place as the longitudinal component of 
the spin current. After taking away the longitudinal com- 
ponent, the remaining part is transverse and therefore sat- 
isfies the Ward-Takahashi identity. The current -magnon 
two-point function which appears in the second term con- 
sists of fermionic and bosonic loop. The latter contains 
magnons and the Higgs boson. 

The Goldstone bosons (f>i and (p2 correspond to the 
SU(2) rotation perpendicular to the local orientation of 
spin (which is along z). 



exp 



(40) 



and they correspond physically to the magnons. / is the 
Goldstone boson form factor which, by virtue of eqn. 
is calculated as the strength of the current-Goldstone- 
boson two-point amplitude. 



3.1 The two-point function and the coupling with 
fermions 

We have noted in eqn. (O that there is a residual symme- 
try associated with U(l)^, which is conserved. The states 



can be classified according to the charges under this rota- 
tion group. First, we define the charge of V't to be +1/2. 
The remaining charges follow automatically, and wc ob- 
tain the values listed in tab. |31 These are necessarily con- 
served. Note that the U(1)em charges are e for the elec- 
tron/hole fields and for all others. These charges are also 
conserved. 



field 


U(l)z charge 


U(1)em charge 


^+ = ip\ 


+ 1/2 


e 




-1/2 


e 


<j)+ = + 102 


+ 1 





<!>- = <I>1 + i(j)2 


-1 





ho = h 









Table 3. The U(l)z and U(1)em charges of the fields, e is 
positive for holes and negative for electrons. 



In order that the Ward-Takahashi identity is satisfied 
by the vertex of eqn. (15^ . the following identity needs to 
be satisfied: 



xqf" 



-fDl\q). 



(41) 



Here, / is a constant of proportionality, and is the same 
quantity as that which appears in eqn. (1^0)) . q^ is the 
momentum fiowing into the two-point function from the 
current (i.e., left to right). This present definition of / is 
more rigorous. Note that this also fixes the sign of D^. Our 
present definition corresponds to taking the couplings to 
be real and taking the sign of D^p to be opposite to that 
for scalar particles. 

Given this definition of /, we can determine the cou- 
pling constants with the fermions by the condition that 
eqn. (|39p satisfies the Ward-Takahashi identity. We then 
obtain the Feynman rules that are given in figs. [2^ and b. 





Fig. 2. The Feynman rules for the coupling of the magnons 
and the Higgs boson with the fermions. 
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For example, for the configuration of fig. [5^, eqn. p9l) 
yields the following Ward-Takahashi identity: 

q^r^ + i\-fD^\q))D^{q){f-'AE)=Gz'-G+\ (42) 

Note that the same Feynman rules can be obtained, less 
rigorously, by considering the rotation associated with the 
Goldstone bosons, eqn. (HO]), and considering its coupling 
with the fermions in eqn. 

The vertices that involve the Higgs boson, which are 
shown in fig. [2j; and d, cannot be fixed by this particular 
type of Ward-Takahashi identity. However, they can be 
fixed by considering the current insertion in the three- 
point amplitude, for example, as shown in fig. [3] 





-1/2 



Fig. 3. The three diagrams whose sum must satisfy the Ward- 
Takahashi identity. The crosses correspond to the modified cur- 
rent vertex defined by eqn. ()39|) . 



The Ward-Takahashi identity applied to fig. [3] also al- 
lows us to determine the magnon-magnon-Higgs vertex. 
However, for doing so, we need to know the form of D^{q) 
and Dfiiq). Let us therefore calculate the current -magnon 
two-point function of eqn. (HH. 




f-'AE 



(b) 



+1 



-T/2 
Fig. 4. The fermionic (a) and bosonic (b) contributions to 
the current -magnon two-point function. 



For now, we calculate the fermionic loop, which is 
shown in fig. 2^. It should be noted that the end result 
of this calculation is independent of whether we consider 
0+ or Using the Feynman rule of fig. [5^, we obtain 



^•^two-point(9) 



In particular, 
expression: 



for the case g — 0, we obtain the exact 



/rjd+ll. 
^---^j-^r^[G^{k)-G+{k)], 

(44) 

where we made use of eqn. p7D . For Fermi liquids, = 1 
and, by symmetry, the spatial components of this ampli- 
tude usually vanishes at g 0. We thus obtain 



4^ 



two— point 



i„t('Z^O) = ri(p_-p+, 0,0,0). 



(45) 



Note that this vanishes for the case of anti-ferromagnetism 
where there is no global spin asymmetry. 

When we compare eqn. (145p with eqn. (|^T|). we see 
immediately that 



~ — p+ (ferromagnetism) , (46) 



and 



-D^^{q) = qo — /2m^ (ferromagnetism), 
-D'^^{q) = — u^q^ (anti-ferromagnetism), 

for small energy and momenta. The inclusion of the bosonic 
loop does not alter this conclusion. The unusual negative 
sign of D refiects the fact that the magnons are pseudo- 
scalar. That is, the fields are rather than (f> in our con- 
vention. 

We need to calculate by other means, such as cal- 
culating the two-point amplitude for finite q, for anti- 
ferromagnetism. u and are parameters which are in 
principle calculable by, for example, evaluating the finite 
q case. However, the form of eqn. (j43p implies rrie 
and u ^ vp. 

There is actually a smarter method than to calculate 
the finite-(7 case (which is cumbersome), but the full cal- 
culation, in the case of anti-ferromagnetism, requires our 
knowledge of the bosonic three-point functions. Let us 
therefore postpone the calculation of anti-ferromagnetic 

and other parameters for now. 

Let us summarize the results of this section up to here. 

Firstly, we summarize the propagators and the two- 
point functions in fig. [S] The Higgs-boson Green's func- 
tions are given in fig. [5}), with a constant energy gap Ah- 
Ah is defined as the energy gap for ferromagnetism and the 
energy gap squared for anti-ferromagnetism. This defini- 
tion is convenient when we discuss the bosonic three-point 
functions. Note that it is an approximation to say that 
Ah is independent of momenta and energy. However, it 
becomes easier to implement current conservation in this 
manner. For completeness's sake, we also list the screened 
photon Green's function (with the approximation that the 
screening, TTem, is constant) and the Higgs-photon mix- 
ing. These are as given in the previous section. 

Secondly, the fermionic vertices are as given before in 
fig. [21 We did not list the photonic vertex, but this is 
given by e. Note that the couplings given in fig. [2] can be 
summarized in the following compact form (c.f. eqn. (j34p '): 



4ff = (/-Mi?)V'^<?.a^, 



(48) 
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(a) 
<i 



F- ^ 

?0-q /2mj,+i0 



AF 



(c) 

X - " - 4>± 

F: (l,q/m^)/ 

AF: {qo,u\)} 

(d) 
1 



+1 



(b) 



F: 



AF: 



95- 



(e) 

7 ^\/\/»- 

F: -e(/-MS)7I,^ixin 



Fig. 5. The propagators and the two-point functions. 77em 
is given by the — component of eqn. (|13|) . and TTmixing is 
given by the — component of eqn. ()10|) . 



and (73 is for the —1/2 fermion. We denote the momentum 
which goes into the vertex by q, so that 9 + 91 = (?2 + 93- 
It is not necessary that the fermions and the bosons are 
on shell, i.e., G^^(gi) etc. need not be zero. 

Upon contraction with g, the first two diagrams yield 



g^^^ = -i^r'AE (1 - G;i(gi)G_(g + qi)) (53) 



and 



(54) 



The amplitude as a whole satisfies the Ward identity if 
the third amplitude satisfies 



q^A^, = 2l^r^AE (1 + Dy^\q2)D^{q2 - q)) . 
This requires vertices of the form shown in fig. [S) 



(55) 



(a) 



(b) 



X 



±1 



±1 / ^ ±1 



qi r q2 



-2(1 



qi+aa. ' 



2m ' 

AF: -2(gi +92) 



Fig. 6. The bosonic three-point functions. In (b), gi -I- 52 is a 
short-hand notation for ({qi + g2)o,'U^(qi + q2)). 



where is defined by 

<P^{^i,h,-v + h). (49) 
(f>i and (f)2 are as shown in tab. [31 and are given by 



62 = -(0++0-). (50) 



u is a parameter, which has the interpretation as the vac- 
uum expectation value of the <P field. In order that the 
energy difference between the two states that is given by 
eqn. (|48p should agree with the actual energy difference 
AE, V needs to satisfy 



V = //2. 



(51) 



^ is essentially a magnetic order-parameter field. This dif- 
fers from the more conventional form such as 



ui(f>i,4>2){o,o,v + h)'' 



(52) 



but they match in the limit of small fields, up to some 
differences in convention. 



3.2 Bosonic vertices 

Let us consider the Ward-Takahashi identity correspond- 
ing to the amplitude described by fig. [3l 

We denote the initial state momentum to be qi and the 
final state momenta to be q2 and q^ . 92 is for the ho boson 



(a) 

f 

+1 ■ 



(b) 

-1 X " ^0 



r 



+1 



X 



(c) 



+1 
- ^ 



+1 



±1 







(d) 



X 



±1 



Fig. 7. The three diagrams (a-c) whose sum must satisfy the 
Ward-Takahashi identity. The bosonic four-point function (d) 
is fixed as a result. 



Finally, we require the Ward-Takahashi identity for 
the sum of the three diagrams which are shown in fig. [7^- 
c. We choose the momenta to be h{qi) -> h{q2) + (f>+{q3), 
with q = q2 + qa — qi being the four-momentum flowing 
into the current. We obtain 

q^Ai^ = -Ai^f-'AhiD^\q,)D^iqi + q) + 1), (56) 



and 



-Ai^r^Ah{Dl\q2)D^{q2 -q) + 1). (57) 



10 



Kosuke Odagiri: Current conservation and ratio rules in magnetic metals with Coulomb repulsion 



Thus we require 



(58) 



This is divergent at large momenta, and needs to be cut-off 
at |k| = K, where K ^ n/a. We then obtain 



in order that the Ward-Takahashi identity is satisfied. 
Hence we obtain the Feynman rule shown in fig. [7Ji. 



4 Calculation of the parameters 

In the preceding section, we worked out the form of the 
theory. Let us now work out the parameters. 

In sec. [21 we used the condition of tadpole cancellation 
to work out a certain rule involving the ratios of density 
of states, that need to be satisfied in ferromagnetism. 

In sec. 121 we presented the Feynman rules and were 
able to relate the current-magnon two-point function to 
the form of and the bosonic propagators. 

We now generalize these results, and work out the four 
parameters, which are (1) (2) AE, (3) Ah and (4) u 
or m^. 

Corresponding to these four unknowns, we have four 
equations, which involve: (1) tadpole cancellation, (2) the 
time component of the current-boson two-point function, 
(3) the space component of the same two-point function 
and (4) the Higgs-boson self-energy. 



4.1 Tadpole cancellation 

Let us start with the condition of tadpole cancellation. 

We treated the ferromagnetic case in sec. [21 The anti- 
ferromagnetic case is calculated analogously, but the mech- 
anism of cancellation is different. This time, we have the 
contribution of the magnon loop: 



PM \AE\ 
2Ah 



K/Attu (d = 2), 
K'^/^^^u (d = 3). 



(64) 



Note that the propagators are all-order, and this requires 
that AE is bare, and that the vertex corrections are not in- 
cluded in eqn. (j59D . Whether AE is stable against higher- 
order corrections depends on the size of the coupling f~^AE 
and the relative size of u compared with the electron ve- 
locity V vp. The form of eqn. (|5^ corresponds to the 
all-order vertex, and therefore this equation, and eqn. (|64p 
which follows from it, suffer from double counting. How- 
ever, this ambiguity, that is due to double counting, can- 
cels when we discuss the Higgs-boson self-energy later on. 



4.2 Current-magnon two-point function 

In sec. [3l we calculated the fermionic contribution to the 
current-magnon two-point function, which is shown in 
fig. [3^. We obtained 



fermionic 



(27r)° 



'^'^^{r'^E)r^G+{k)GAk- q). 

(65) 

For the simple case of e(k) — (fi,k)^/2TOe, the vertex is 
given by — (1, (k— q/2)/me). The bosonic loop, which 
corresponds to fig. [Jb, is written as 



bosonic 



^^p^(2/-iZ\,0(-2(2fc - q))Du{k)D^{k - 



^^matnon = ^'(2/-^,) 



which must be equal and opposite to the fermionic loop: 

^|^(G+(fc)-G_(fc)) 
. . 

Note that the sum over positive and negative AE is im- 
plicit. 

This gives us the following condition: 



(27r) 



(61) 

Note that AE is positive when p_ > p+. That is, the right- 
hand side is negative. Let us introduce a more compact 
notation: 



^ ^ (j^) (KQ) This is for the anti-ferromagnetic case. Note that the case 

{2tt)'^+^ of ferromagnetism requires a further twist, as we have not 

yet included the higgs-screened-photon mixing. 

Let us consider the limit of small external momentum, 
q — > 0. It is easy to see that the fermionic amplitude van- 
ishes for anti-ferromagnetism. As for the bosonic ampli- 
tude, this vanishes for ferromagnetism because of the ab- 
sence of negative energy states. The amplitude vanishes 
for anti-ferromagnetism also, but for a different reason, 
namely symmetry. 

Let us, instead of trying to evaluate these integrals for 
arbitrary values of q, make use of the Ward-Takahashi 
identities to replace the current-magnon two-point func- 
tions with the corresponding current-current two-point 
functions (c.f., ref. [5]). 

To do so, we first write down the current-current two- 
point functions as 



sublatticcs 



2 Ah J D^^pm\AE\ 



(62) 



fermionic 



The convention is that pM = IP- + P+l is positive. The 
integral of is evaluated easily: 



and 



Da = 



d'^+^k 



-1 



d'^k 



(27r)'i+ii fc2 _ ^2^2 + io J 2it(27r)'« |k| ' 

(63) 



j^^^r^r'^G^G^ik-q), (67) 



(27^)'^^ 



(-2(2fc - g)^)(-2(g - 2kY)Dh{k)D^{k - q). (68) 
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By virtue of the Ward-Takahashi identities, we obtain 

q^ll J y\ — ^formionic 

where C are given by 



'^bosonic (69) 



fcrmionic 



and 



^01^^r^{k,q-k) (G+(fc) - G_(fc - q)) 

(70) 



bosonic 



^^-^^{-2{2k-qr) {Duik) + D^{k q)) . 

(71) 

We now equate A with the Feynman rules of fig. [St, 
and take the derivative with respect to q\ in the limit of 
small q. For the ferromagnetic case, we obtain 



fern 



/2TO^Miag(0,-/) = T^Cf'^rmionic 



d_ 

dqx 



q^O 

(72) 

/ stands for the spatial identity matrix. The — com- 
ponent of this equation is zero on the right-hand side and 
in the second term of the left-hand side, whereas the first 
term on the left-hand side is non-zero: 



00 
fern 



G-i-G ^ 



p- 



AE 



(73) 



This happens because of the approximation ^ [q) — 
go — I'i'rn^. There is, in principle, a q^ term as well, 
the omission of which is inconsistent with the — com- 
ponent of this equation. As for the spatial components, we 
obtain 



{G+{k) 



d'^+^k 
G-{k)) 



2top 



G+{k)G^{k)vl 



(74) 



Here, d/nie refers to the second derivative of e(k), whereas 
Ve refers to the first derivative. 

Let us introduce the following shorthand notation: 



.2 -1 / P- + P+ 



jp- - p+)K 

AEd 



(75) 



There is an obvious generalization to the case of spatial 
asymmetry. We expect this final result to be stable against 
higher-order corrections, because the renormalization fac- 
tors due to the vertex correction and the Green's functions 
cancel. 

As discussed in sec. [51 (p_ + p+) is not well-defined. 
However, the ratio of (p_ + p+) against (p_ — p+) is well 
defined because of eqn. Furthermore, p is given by 
P- - P+- 

As an order estimation, we can say that Ve can be taken 
to be almost constant near the Fermi surfaces. It would 
be a bad approximation to say that mg is also constant, 



but we can introduce a quantity rrie to be the inverse of 
the average inverse fermion mass. We then obtain 



2me 1 2meVp 



P 



AEd 



(76) 



P is the spin asymmetry, is necessarily positive, but 
needs not be positive although we generally expect it 
to be. If TTie is positive, then the inequality reads 



P < 



AEd 

2meVp 



(77) 



Let us now turn to the anti-ferromagnetic case. Here 
we need both the fcrmionic and the bosonic loops. Corre- 
sponding to eqn. (1721) . we now have 



(78) 



The fcrmionic contributions are as given above. The bosonic 
contributions arc given by 



71, 



bosonic 



-16 Jk^k^DnD^, 



and 



d 
dqx 



Cb^osonic = 2diag(l, ~u'l) j Dh~D^ 



(79) 



(80) 



using the same notation as in eqn. (|62[) . The integral over 
is given by eqn. (|63)) . and is a positive quantity. The 
integral over Dh is given by 



d'^'k 



2{2TTfy/u^\i^ + Ah 



(81) 



and this is a negative quantity. 
Altogether, we obtain 



and 



2 j {Dh - D^) -16 J klDhD^ 



Pm 
\AE\ 



(82) 



/2 + 2 y" {Dh - ^0) - ^ / ^^DhD^ 



LIP 

-7/2 



2to„ 



\AE\d 



(83) 



4.3 The Higgs-boson self-energy 

We now come to the final condition, namely that the 
Higgs-boson excitation energy A^ is given by the self- 
energy diagrams which are shown in fig. [5) 

The fcrmionic contribution is similar to 77^" which was 
calculated in sec. [5J and is given by 



(G+(i-)G+(t - 5) + G_(ii)G_(f: - 9)) . (84) 
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(a) 



(b) 



(c) 



Fig. 8. The diagrams for the self-energy of the Higgs boson. 



Hence 



n^ = {r'AEfi9-i^F) + g+{eF)) 



(85) 



Note that AE refers to the bare quantity. This is because, 
firstly, {AE)'^ in eqn. needs to be the product of the 
bare AE and the renormalized AE. However, the renor- 
malization of AE gives rise to the renormalization factor 
which is opposite to the renormalization factor Z for 
each propagator. It follows, therefore, that AE in eqn. (j85l) 
actually refers to the bare quantity. 

At g = (and at zero temperature), the contributions 
of fig. [8)3 and c are zero for ferromagnetism. Hence, for 
the case of ferromagnetism, we obtain 



A,= '-M±9±M^AEr, 

P- - P+ 



(86) 



If the density of states 5 is a linear function, then Ah = 
2AE since (p_ — p+) is given by the area of a trapezium 
whose two parallel sides are g- and g+, and whose height 
is AE. If not, and g is a convex function in between g- and 
g+ as is the case for iron [TT], Ah will be less than 2AE. 
This gives a useful estimate of the Higgs-boson excitation 
energy, which can be tested experimentally. 

We should remember that the Higgs-screened-photon 
mixing cannot be neglected when the spin asymmetry P is 
large. The actual value of Ah where the resonance occurs 
will be sensitive to the behaviour of the photonic modes 
(screened photon and plasmon). 

Let us now turn to anti-ferromagnetism. The contri- 
bution of fig. \Ejp is given by 



- 2i7^ = 



(27r)° 



ji2r'AhrD4k)D^{k-q). (87) 



This is divergent, but is imaginary at zero temperature for 
Qq — u^q^ > (which is where the Higgs mode needs to 
exist). We therefore omit this contribution for now. The 
contribution of fig. [5t is given by 



in!} = 



-jd+ll 



(27r)" 



-{-4f-'Ah)D^ik). 



Now using eqn. ((62|) . this reduces to 

n^ = -Ar^Ahj D^{k) = -2r 



\AE\ 



(89) 



Hence, 

Ah - n1 + n'^ 

= 2f-^ \AE\ 



{g^ieF)+g+{e))\AE\-p,, 



We thus obtain the anti-ferromagnetic ratio rule: 
g-{€F) + g+iep) 



2pm/ \AE\ 



> I. 



(90) 



(91) 



This is satisfied if the density of states is a concave func- 
tion in between the two Fermi energies. 

Let us denote the concavity by Sc, defined as: 



g-{<^F) + g+iep) 



2pm/ \AE\ 



1. 



This then leads to 



2pM \AE\ 



By eqn. (j64L we then obtain 



5cK/nu (d = 2), 
ScK^/2n^u {d = 3). 



(92) 



(93) 



(94) 



Small 6c therefore leads to strong coupling f^^AE. We 
expect physically that strong coupling tends to suppress 
magnetism, because the oscillations between the two spin 
states will become more frequent. Our results are not af- 
fected so long as the densities of states can be defined. 
However, AE will receive a large correction through the 
electron self-energy. 



4.4 Summary of results at zero temperature 

Let us summarize our results. 

For the case of ferromagnetism, the parameters AE, 



/ , and Ah are fixed by the following constraints: 

p+/p- = g+{f-F)/g-{£F), 
f = P-- p+, 

f2 -I _ / P- + P+ jP- - P+)vl 

"""^ " \ 2me AEd 
fAh = {g-{ep)+g+{ep)){AEf. 



(95) 
(96) 

(97) 

(98) 



Out of these equations, which are all non-perturbative, the 
first three are relations that only involve all-order quanti- 
ties. In the last equation, AE refers to the bare quantity. 
In eqn. (I97p . AE is the all-order quantity, but is assumed 
to be more or less independent of energy and momenta 
(though generalization is possible). 

For the case of anti-ferromagnetism, \AE\, u and 



Ah are fixed by eqns. (gll), (|H31) and (1131) ■ Equations 
([M)) and ([M)) involve \AE\ as a bare quantity and Ah 
in eqn. is ambiguous. Equations and only 
involve all-order quantities, but are dependent on the UV 
cut-off, as is the case in eqn. (|64|. 

We obtained the rule (5c > 0, where 5c is a measure of 
concavity and is defined by eqn. (j92p . 
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4.5 Finite temperature analysis 

Since our results involve diagrams that are evaluated for 
g = 0, it is, in principle, straightforward to generalize 
them to finite temperatures. However, the bosonic dia- 
grams, which were zero in the case of ferromagnetism, be- 
come non-zero at finite temperatures, and therefore the 
resulting expressions are messy. 

A full calculation is beyond the scope of this present 
analysis, but let us present two representative results. 

First, for the case of anti-ferromagnetism, we have 
found that the bosonic loop of fig. [Hd is real and diverges 
for finite T: 

. ^ {2r^Ahf f ^ ^ / coth(x) \ 

This makes Ah negative, and so magnetic order is forbid- 
den. In our opinion, this implies that in anti-ferromagnetic 
metals, a genuine long-range order is not permitted, at 
least at finite temperatures. 

Second, let us consider how the ratio rule of eqn. ([T^ is 
modified at finite temperatures. We now have the bosonic 
contribution which reads 

__!^^ t 

(27r)'i exp((k2/2m^)/T)- f 

(100) 

This is evaluated easily using standard methods. For the 
case of three spatial dimensions, we obtain 

^^oltiT) ^ -2/-iZi,C(3/2) (^)"' ■ (101) 

Here, C(3/2) = 2.612 ••• . This contribution should be 
equal and opposite to the fermionic contributions, which 
are given by 

^tr^liT) = ir'AE) [-{p, p^)t + {p, + Pt)TP(T)] . 

(102) 

Here p and P correspond to their finite-temperature coun- 
terparts: 

PT = j f{{e-p)/T)g(e)de, (103) 
gr^- I riie-fi)/T)g{e)de, (104) 
where / is the Fermi distribution function. Hence 

- {Pi - Pt)t + {Pi + PthPiT) = 2C(3/2) 

(105) 

Ah, and AE are also functions of temperature. 

For small T, we can assume that only the first term 
on the left-hand side depends significantly on T, and that 
the parameters on the right-hand side can be taken as 
constants. This then implies that the magnetization goes 
down as T^/^, which is a well-known result. All of the 
parameters on the right-hand side are, in principle, mea- 
surable. This can then be tested experimentally. 



5 Conclusions and Outlook 

We presented a nonperturbative framework for treating 
magnetic order in metals, caused by a Coulomb interaction 
(or generalized Coulomb interaction). 

We obtained interesting 'ratio rules' involving the den- 
sities of states for both ferromagnetic and anti- ferromagnetic 
cases. These involve all-order quantities (with the excep- 
tion of l^i^l) and can therefore be compared directly with 
the experimental numbers, if they become available at 
greater precision. 

We have seen that the shape of the density-of-states 
curve play an essential role in determining the possibil- 
ity of magnetic ordering. The density of states must rise 
■with energy, when the charge carriers are electrons, for 
ferromagnetism. 

For anti-ferromagnetism, the density-of-states curve must 
be concave. However, we have seen that the radiative cor- 
rections, due to the magnons, at finite temperatures breaks 
long-range orders. In our understanding, this means that 
genuine long-range anti-ferromagnetic order is not possi- 
ble, at least at finite temperatures. More work is required 
to elucidate the nature of the ground state. 

The case of magnetic insulators is not covered by this 
work, in which the exchange energy AE is considered to 
be more or less independent of k. 

Two cases require special attention, which we have not 
been able to discuss in much detail. The first is the case 
of strong coupling, which occurs when /, or the vacuum- 
expectation value V of the magnetic order-parameter field, 
is small. Here, we expect that the radiative corrections 
suppress the magnetic order and that the system will favour 
the paramagnetic state. The second is the case of large 
magnon velocity u, in comparison with the electron ve- 
locity wf, in the case of anti-ferromagnetism. Here, the 
response of the magnetic background becomes instanta- 
neous towards the movement of the electron. We hope to 
be able to discuss this case in a separate publication [T3] 

The results of this work can be used to calculate arbi- 
trary amplitudes, such as scattering amplitudes. 

The methods presented in this work, being an adap- 
tation of Gribov's analysis of axial-current conservation, 
is of a general nature. However, we are presently unaware 
of other possible applications of the methods presented 
herein. 

We thank I. Hase, S. Sharma, K. Yamaji and T. Yanagisawa 
for extensive and informative comments and discussions. 

We have been informed by Dr. I. Hase that a phenomeno- 
logical study of the correlation between densities of states of 
a material and its magnetic properties has previously been re- 
ported. However, we have not been able to locate this study. 
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